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Abstract — A distributed MIMO radar is considered, in which 
the transmit and receive antennas belong to nodes of a small scale 
wireless network. The transmit waveforms could be uncorrelated, 
or correlated in order to achieve a desirable beampattern. The 
concept of compressive sampling is employed at the receive 
nodes in order to perform direction of arrival (DOA) estimation. 
According to the theory of compressive sampling, a signal that 
is sparse in some domain can be recovered based on far fewer 
samples than required by the Nyquist sampling theorem. The 
DOAs of targets form a sparse vector in the angle space, and 
therefore, compressive sampling can be applied for DOA esti- 
mation. The proposed approach achieves the superior resolution 
of MIMO radar with far fewer samples than other approaches. 
This is particularly useful in a distributed scenario, in which the 
results at each receive node need to be transmitted to a fusion 
center. 

Keywords: Compressive sampling, MIMO Radar, DOA Esti- 
mation 

I. Introduction 

A multiple-input multiple-output (MIMO) radar system, as 
originally proposed in [l]-[5] transmits multiple independent 
waveforms via its antennas. Consider a MIMO radar equipped 
with Mt transmit and M r receive antennas that are close 
to each other relative to the target. The phase differences 
induced by transmit and receive antennas can be exploited 
to form a long virtual array with M t M r elements. This 
enables the MIMO radar system to achieve superior spatial 
resolution as compared to a traditional radar system. MIMO 
radar transmitting correlated signal waveforms in order to 
achieve a desired beampattern has also been proposed [10]- 
[12]. This is useful in cases where the radar system wishes 
to avoid certain directions, because they either correspond to 
eavesdroppers, or are known to be of no interest. 

Compressive sensing (CS) has received considerable atten- 
tion recently, and has been applied successfully in diverse 
fields, e.g., image processing [6] and wireless communications 
[7]. The theory of CS states that a A'-sparse signal x of 
length N can be recovered exactly with high probability from 
(D(K\ogN) measurements via linear programming. Let ^ 
denote the basis matrix that spans this sparse space, and let 
$ denote a measurement matrix. The convex optimization 
problem arising from CS is formulated as follows: 



|s||i, subject to y = $x = <E>^I/s 



(1) 



where s is a sparse vector with K principal elements and the 
remaining elements can be ignored; $ is an M x N matrix 
incoherent with ^ and M <C N. 



In this paper, we propose a distributed MIMO radar system, 
where transmit and receive antennas belong to nodes of a 
wireless network that are uniformly distributed on a disk or a 
certain radius, located wireless network nodes. The readings 
of the receive nodes are transmitted to a central node for 
DOA estimation. Energy efficiency is an important issue in 
such a wireless network as the nodes operate on battery. We 
employ the idea of compressive sampling in order to save in 
energy consumed during data transmission to the central node. 
Recently, the work of [8] considered DOA estimation of signal 
sources using CS. In [8], the basis matrix ^ is formed by 
the discretization of the angle space. The source signals were 
assumed to be unknown, and an approximate version of the 
basis matrix was obtained based on the signal received by a 
reference vector. The signal at the reference sensor would have 
to be sampled at a very high rate in order to construct a good 
basis matrix. Here, we extend the idea of [8] to the problem of 
DOA estimation for MIMO radar. Since the number of targets 
is typically smaller than the number of snapshots that can be 
obtained, DOA estimation can be formulated as the recovery 
of a sparse vector using CS. Unlike the scenario considered 
in [8], in MIMO radar the transmitted waveforms are known 
at each receive antennas. This enables each receive antenna to 
construct the basis matrix locally, without knowledge of the 
received signal at a reference sensor or any other antenna. We 
consider the more general case of correlated signal waveforms. 
We provide analytical expressions for the average signal-to- 
jammer ratio (SJR) for the proposed approach. Simulation 
results show that the proposed approach can accomplish the 
super-resolution of MIMO radar systems while using far fewer 
samples than existing methods, such as Capon, amplitude and 
phase estimation (APES) and generalized likelihood ratio test 
(GLRT) [2]. In particular, the proposed approach can enable 
each node to obtain a good DOA estimate independently. 
Further, it results in much less information to be transmitted to 
a fusion center, thus enabling savings in terms of transmission 
energy. 



II. Signal Model for MIMO Radar 

We consider a MIMO radar system with Mt transmit nodes 
and M r receive nodes that are uniformly distributed on a disk 
of radius r. For simplicity, we assume that targets and nodes 
lie on the same plane. Further, we assume that each node in 
the network knows which are the nodes that serve are transmit 
and receive antennas and what their coordinates are relative to 
a fixed point in the network. This information can be provided 
by a higher network layer. Let us denote the locations in 
rectangular coordinates of the i-th transmit and receive antenna 
by (x\,y\) and (x\,y\), respectively (see Fig{T]>. 

The location of the fc-th target is denoted by the po- 
lar coordinates (dk,9k), where dk is the distance between 
this target and the origin, and 9 k is the azimuthal angle, 
which is the unknown parameter to be estimated in this 
paper. Under the far-field assumption dk S> J {x\) 2 + {y\) 2 
and dk -J {x\) 2 + (y[) 2 , the distance between the ith 
transmit/receive antenna and the fc-th target d\ k ld r ik can be 



d k - v/ r (h), where r/f r (0 fc ) 



approximated as d l J k r 
xl /r cos(6> fc ) + y\ lr sin(0 fc ). 

Let Xi (n) denote the discrete-time waveform transmitted by 
the i-th transmit antenna. Assuming the transmitted waveforms 
are narrowband and the propagation is non-dispersive, the 
received baseband signal at the fc-th target equals [4] 



Vk{n) 



Pk^2xj(n)e 

i=i 



= l3 k e-^ d ^ T {n)w{9 k ) k = l,...,K (2) 
where A is the transmitted signal wavelength, 

v(0 fc ) = [e^^\...,e j ^^} T (3) 
x(n) = ...,x Mt (n)] T . (4) 

Due to reflection by the target, the Z-tii antenna element 
receives 

K 

zi(n) = ^e-^S(»)+6,(n),l = l,..,M r (5) 

fe=i 

where ei (n) represents independent and identically distributed 
(i.i.d.) Gaussian noise with variance a 2 . 

On letting L denote the number of snapshots, we have 

*l(0) 



The ft-th target 



The /-th transmit 
anntenna / 
The i-th receive / / 
anntenna / 



where 



Fig. 1. MIMO Radar System 
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e~°~ k /3k if there is target at a n 
otherwise 



III. Compressive Sensing for MIMO Radar 

Assuming that there exists a small number of targets, the 
DOAs are sparse in the angle space, i.e., s = [si, . . . , sjv] is a 
sparse vector. A non-zero element with index j in s indicates 
that there is a target at the angle a,. 

By CS theory, we can construct a basis matrix for the 
Z-th antenna as 



*i = [e J T"r(«i)Xv( ai ), . . . , e^ r(QN) Xv(aAr)] 



(8) 



Ignoring the noise, we have z; = ^;s. Then we measure 
linear projections of the received signal at the Z-th antenna as 



(9) 



where is an M x L random Gaussian matrix which has 
small correlation with Combining the output of N r receive 
antennas, we have 



z/ 



zi(L-l) 



K 


ri 




$1*1 










fc=l 













s, 1 < Nr < M r 



(10) 



K 



J2 e-W d *- T M°>»0 k X.v(9 k ) + e, 



(6) 



fc=i 



where y fe = [y k {0), y k {L - 1)] T , e z = [e,(0), . . . , e,(i - 
1)] T and X= [x(0),...,x(i- 1)] T . 

By discretizing the angle space as a = \ot\, . . . ,ajv], we 
can rewrite (0 as 



Therefore, we can recover s by applying the Dantzig selector 
to the convex problem in ( [Tol l as in [9]: 



s = min ||s||i s.t.\\<d H (r — 9s)||oo < fx. 



(ID 



Z[= ^ e3 f,[K) SnXv(an)+e! 



(7) 



71=1 



According to [9], we can recover the sparse vector s with very 
high probability if we select fi = ( 1 + 1 ~ 1 ) y2 log Na 2 , where 
t is a positive scalar and a 2 is the noise power. 



IV. Performance Analysis in the presence of a 

JAMMER SIGNAL 

In the presence of a jammer at location (d, 9) the signal 
received at the Z-th receive antenna can be represented as 

K 

n = $ ; ]T e -^ 2rffc -^»/3 fc Xv(0 fe ) 



fc=l 



+$ ie - J ^ (d ~"' r(e)) /3b + $, 



(12) 



where (3, b denote respectively the reflection amplitude and 
waveform of this jammer. Since b is uncorrelated with the 
transmitted waveforms X, the effect of the jammer signal is 
similar to that of addictive noise. Let A; = and D; = 

X ff A;X, where D(i,j) denotes the (i,j)th element of D. 
We assume that the TX/RX nodes are uniformly distributed 
on a disk with the radius r. Thus, the average power of the 
desirable signal P s (l) can be represented by 

K 



P S {1) = E{ V e j ^ [2{dk ~ dk ' ] ~ W ( 9k > ~ r >< < 9 *' > \ (3* k (3 k , 

k,k' = l pi(kM') 

K 

x v H (9 k )X H A l Xv(9 k ,)} = E{Y] \(3k\ 2 Qkk} 

+E{J2 Pi(Kk')PlPk>Qkk>} (13) 

where Q w = . D,(i, j^W^-vt^)). 

Following [3], we know that fh(h) = — V^~- ~ h?, — 1 < 
h < 1 if h = jsm($), where f f (f) = 2j,0 < f < r and 
/^(^) = _7T < V < 7T. Then ^-jv^} = 2^*, where 
Ji(-) is the first-order Bessel function of the first kind. Using 
this property and letting ay = r] t / r (9 k i) — vt^i^k), we have 



£{e^ 0i O = l ^HV 1 )) *=3 ( i4) 
1 1 I rf(2) i^j (i4) 

where t^(.t) = 2 Jl ^ r A . 

Therefore, the average power of the desirable signal P s (l) 
taken over the positions of TX/RX nodes can be found to be: 

K 

P S {1) = £ |ft| 2 ED 1 (M)+^D i ( lJ V(2)] 



fc=i 



£ /3^e^^-^)r,(4sin(^-^)) 



x [7 7fefe ,yD ; (»,i) + yD i (i,j)?7 2 (2)]. (15) 

Similarly, the power of the jammer signal is given by 

Pj(l) = ( e -i s r(<l-nr(0))p^ e - j ^(d- v r(8)) /3 y 

xb H A ; b = |/3| 2 b H A ; b . (16) 

The SJR given the node locations is the ratio of the power of 
the signal over the power of the jammer. Since the denominator 



does not depend on node locations, the average SJR equals the 
ratio of ( [13] ) and ( fl6b . 

Since the jammer signal is uncorrelated with the transmitted 
signal, the SJR can be improved by correlating the jammer 
signal with the transmitted signal. Combining this with CS, 
the measurement matrix in (O is modified as 



H 



(17) 



Moreover, since $/ is a Gaussian random matrix, <P; is still 
Gaussian; therefore it satisfies the restricted isometry property 
(RIP) and is incoherent with ^>i, thus guaranteeing a stable 
solution to (fTTT i. Based on ( TPTI i. the average power of the desir- 
able signal P s (l) is equal to O except D ; = X H XA ; X H X. 
The average power of the jammer signal using $/ is rewritten 
as Pj(l) = |/3| 2 b ff XA;X ff b. 

Approximating X H X ~ Ij\/ t and using b H b = 1, the SJRs 
based on $/ and $/ can be approximated as — - ffi an d 

l^p , respectively. Therefore, the SJR using ( TT7l > can 
be generally improved by a factor of L/M t since L ^> M t . 
the DOA estimates can be improved by the increase in L. 
However, the time duration of the radar pulse might need to 
be longer as well. 

As simulation results show (see Section [V]), the proposed 
method can yield good performance even using a single re- 
ceive antenna. With a good initial estimate of DOA, the receive 
nodes can adaptively refine their estimates by constructing a 
higher resolution basis matrix ^/ around that DOA. Restricting 
the candidate angle space, may reduce the samples in the angle 
space that are required for constructing the basis matrix, thus 
reducing the complexity of the l\ minimization step. On the 
other hand, the resolution of target detection can be improved 
by taking the denser samples of the angle space around the 
intimal DOA estimate. Furthermore, the transmit node can 
design the correlated waveforms for transmit beamforming as 
well based on the good initial estimate. 

V. Simulation Results 

In this section, we consider a MIMO radar system with the 
transmit/receive antennas uniformly distributed on a disk of 
radius 10m. The number of transmit nodes is fixed at M t = 50. 
The carrier frequency is 8.62 GHz. A maximum of L = 512 
snapshots are considered at the receive node. The received 
signal is corrupted by zero mean Gaussian noise. The SNR is 
set to 20 dB . 

There are two targets located at 9^ = —1°, 1°, with 
reflection coefficients /3k = 1, k = 1,2. A jammer is located 
at 15° and transmits an unknown Gaussian random waveform 
and with amplitude 20, i.e., 26 dB above the target reflection 
coefficients j3k- We sample the angle space by increments of 
0.5° from -8° to 8°, i.e., a = [-8°, -7.5°, . . . , 7.5°, 8°]. 

First, we compare the performance of DOA estimation using 
the proposed method and three approaches [2], i.e., the Capon, 
APES and GLRT techniques. Fig. [2] and Fig. [3] show the 
modulii of the estimated reflection coefficients (3k, as functions 
of the azimuthal angle for N r = 1 and 10 receive antennas, 



respectively. In Fig. [2] we use the uncorrelated QPSK wave- 
forms; while in Fig. [3] we use correlated waveforms designed 
according to the desired beampattern Pd(a n ) as 

PI \-i 1 - 3 °<««< 3 ° nsl 
d[an> _ \ -8° < a n < -3° and 3° < a n < 8° ' U8J 

Based on that beampattern, the method of [10] was followed to 
design R. Then the transmitted waveforms can be constructed 
as x(n) = R5w, where w is a i.i.d random vector with zero 
mean and E{ww H } = I/L. 

In both (a) and (b), the top three curves correspond to the 
azimuthal estimates obtained via Capon, APES and GLRT, 
using 512 snapshots. The bottom curve is the result of the 
proposed approach, obtained using 35 snapshots only. One 
can see that in the case of using only one receive node, 
the presence of the two targets is clearly evident via the 
proposed method based on 35 snapshots only using both inde- 
pendent and correlated waveforms. The other methods produce 
spurious peaks away from the target locations. When the 
measurements of multiple receive nodes are used at a fusion 
center, the proposed approach can yield similar performance 
to the other three methods. However, the comparison methods 
would have to transmit to the fusion center 512 received 
samples each, while in the proposed approach, each node 
would need to transmit 35 samples each. 

The threshold /i in (fTTTt affects DOA estimation for the 
proposed method. The increase in /i while keeping Mt and Nt 
constant can reduce the ripples of DOA estimates at the non- 
target azimuth angles at the expense of the accuracy of the 
target-reflection-coefficient estimates. The increase in /i can 
also reduce the complexity of (TTTT) because the constraint is 
looser than that of smaller /i. If fj, is too large, however, the £\- 
norm minimization does not work. In Fig|2] and [3] relatively 
large thresholds, i.e., /i = 12, 10, were used for the single 
receive node case. As a result, the CS method yielded less 
accurate estimates of the reflection coefficients magnitude than 
the Capon and APES, but with very few ripples. 

Finally, we discuss the effect of L, N t and M t on the 
performance of the Capon, APES, GLRT and CS. Fig. H 
compares the performance of these four approaches using 
independent waveforms for different combinations of N r and 
L, whose product is fixed at 512. In order to quantify the 
performance of DOA estimation, we define the ratio of the 
square amplitude of the DOA estimate at the target azimuth 
angle to the sum of the square amplitude of DOA estimates at 
other angles as the peak-to-ripple ratio (PRR). Fig.[5]compares 
PRR as a function of L for these four approaches using 
uncorrelated signal waveforms. We consider the scenarios in 
which N r = 1,5,10,30. For fixed N t and M t , the increase 
in L can improve the performance of these four methods. In 
the presence of a moderate jammer, APES and CS can yield 
relatively accurate DOA estimates even with a small L. For 
Capon and GLRT, N r must be greater than L in order to obtain 
a nonsingular sample covariance matrix of the received signal 
R. This is because Capon and GLRT need to calculate the 
inverse of R. On the other hand, an increase in both M t and 



N r can also improve the performance of the Capon, APES, 
GLRT and CS while L is fixed. If either M t or N r is too 
small, even a significant increase in the other parameter cannot 
improve performance of the first three approaches. However, 
CS can yield the desired DOA estimates even with a single 
receive antenna with a sufficient L and Nt. In the senecios 
considered in our simulations where L = 2 9 , for instance, M t 
and N r are required to be greater than 8 to yield the desired 
DOA estimates using the Capon, APES and GLRT. With a 
single receive node, CS requires at least 20 transmit nodes, 
while in the cases of multiple receive nodes, the requirement 
of M t and N r for CS is the same as in the other three methods. 

VI. Conclusion 

We proposed a distributed MIMO Radar system imple- 
mented by a small size randomly dispersed wireless network. 
There are several advantages in using the proposed distributed 
approach as opposed to using a standard linear array. The radar 
system can be easily deployed; no pre-existing infrastructure 
is required. In a high density network there are many degrees 
of freedom to design the beampattern as desired around the 
look direction, which is important for clutter reduction or 
for reduction of scanning time. By randomizing the set of 
transmitters and receivers we can used the network power 
efficiently. By selecting well separated nodes we can increase 
spatial diversity. The resolution can be easily adjusted by 
employing more or less transmit nodes. The radar system is 
robust; should some nodes be deactivated the system perfor- 
mance will not be affected. 

For the proposed MIMO radar system, a compressive sens- 
ing method has been exploited to estimate the DOAs of targets 
using both independent and correlated waveforms. The DOA 
of targets can construct a sparse vector in the angle space. 
Therefore, we can solve for this sparse vector by ^;-norm mini- 
mization with many fewer samples than conventional methods, 
i.e. the Capon, APES and GLRT techniques. The proposed 
method is superior to these conventional methods when one 
receive antenna is active. If multiple receive antennas are used, 
the proposed approach can yield similar performance to the 
other three methods, but by using far fewer samples. 
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Fig. 3. DOA estimates of two targets with close azimuthal angles using 1 (left 
column) and 10 (right column) receive antenna using independent waveforms. 
The top three curves were obtained using 512 snapshots. The bottom curve 
was obtained using 35 snapshots only. 
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